Halo nucleus Be: a spectroscopic study via neutron transfer 
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The best examples of halo nuclei, exotic systems with a diffuse nuclear cloud surrounding a tightly- 
bound core, are found in the light, neutron-rich region, where the halo neutrons experience only 
weak binding and a weak, or no, potential barrier. Modern direct reaction measurement techniques 
provide powerful probes of the structure of exotic nuclei. Despite more than four decades of these 
studies on the benchmark one-neutron halo nucleus ^^Be, the spectroscopic factors for the two 
bound states remain poorly constrained. In the present work, the ^°Be(d,p) reaction has been used 
in inverse kinematics at four beam energies to study the structure of ^^Be. The spectroscopic factors 
extracted using the adiabatic model, were found to be consistent across the four measurements, and 
were largely insensitive to the optical potential used. The extracted spectroscopic factor for a 
neutron in a nlj = 2sx/2 state coupled to the ground state of ^"Be is 0.71(5). For the first excited 
state at 0.32 MeV, a spectroscopic factor of 0.62(4) is found for the halo neutron in a lpi/2 state. 

PACS numbers: 25.60.Je, 25.60.Bx, 25. 45. Hi 
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Nuclear halos are a phenomenon associated with cer- 
tain weakly-bound nuclei, in which a tail of dilute nu- 
clear matter is distributed around a tightly bound core 
[IH5] . This effect is only possible for bound states with no 
strong Coulomb or centrifugal barrier, and which lie close 
to a particle-emission threshold. Though excited-state 
halos exist, the number of well-studied halo states is pre- 
dominantly limited to a handful of light, weakly-bound 
nuclei which exhibit the phenomenon in their ground 
state. 

The neutron- rich nucleus -'^^Be is a brilliant example 
of this phenomenon, with halo structures in both of its 
bound states, and light enough to be modeled with an 
ah initio approach. It is well documented that the 1/2+ 
ground state and l/2~ first excited state in ^-'^Be are in- 
verted with respect to level ordering predicted from a 
naive shell model. There has been considerable theoreti- 
cal effort toward reproducing this level inversion in a sys- 
tematic manner, while maintaining the standard ordering 
in the nearby nuclide ^'^C, where the 1/2+ state lies over 
3 MeV above the l/2~ ground state. A Variational Shell 



Model approach |;4 and models which vary the single- 
particle energies via vibrational [15 and rotational [5] 
core couplings reproduce this level inversion in a system- 
atic manner. Common to the success of these models is 
the inclusion of core excitation. Ah initio No-Core Shell 
Model calculations [7] have been unable to reproduce this 
level inversion though a significant drop in the energy of 
the 1/2+ state in ^-'^Be is reported with increasing model 
space. In all of these models, the wave functions for the 
^^Be halo states show a considerable overlap with a va- 
lence neutron coupled to an excited ^°Be(2+) core, in 
addition to the naive n(g)^°Be(0+) component. Despite 
decades of study, the extent of this mixing is not well 
understood, with both structure calculations and the in- 
terpretation of experimental results ranging from a few 
percent to over 50 percent core-excited component. 

Experimentally, it is possible to gain quantitative in- 
sight into the mixed configuration of a state by study- 
ing reactions which provide observables that are sensi- 
tive to different components of the nuclear wave function. 
By comparing the measured differential cross sections to 
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FIG. 1. A selection of literature values for spectroscopic 
factors for the ground and first excited states of ^^Be. 
Experimentally extracted spectroscopic factors are taken 
from references [8](a), l9](b), [l0](c), [II](d,e), [l2](f), [T3](g), 
[14] (h), and [l5](i). Results of theoretical calculations are 
taken from [5](j), [6](k), [7](1), and |4](m). Filled black 
symbols denote spectroscopic factors extracted from (d,p) 
data. 



those calculated theoretically, a spectroscopic factor S 
can be extracted, which reflects in a model-dependent 
way the extent to which the nuclear state studied re- 
sembles that used in the calculation. Spectroscopic fac- 
tors reported from numerous direct-reaction studies, in- 
cluding one- neutron transfer [H [9l [131 IE] > two-neutron 
transfer [T7], neutron knockout [13], Coulomb breakup 
[T5] , and reanalyses of the neutron-transfer data [TUl [H] , 
are shown in Figure [T] along with those from theoretical 
calculations. For the ground state, the spectroscopic fac- 
tors correspond to overlaps with a 2si/2'8)^''Be(0+ ) con- 
figuration. For the first excited state, they refer to the 
lPi/2'8i"'^°Be(0i,). These values exhibit significant varia- 
tion beyond the reported uncertainties, even within the 
same reaction type. 

The challenges associated with determining spectro- 
scopic factors from measurements using various reac- 
tions has received considerable attention in recent years 
[T51 [T^ . Even extracting consistent spectroscopic factors 
from a single reaction type is sensitive to the description 
of the reaction mechanism (evidenced by the variation in 
S for different analyses of {d,p) data), and to the nuclear 
structure included in the model (evidenced by the range 
of S from {p,d) studies [131 US] using different potentials 
to calculate the wave functions). However, it should be 
expected that consistent results are obtained for mea- 
surements employing the same reaction type, provided 
the data are analyzed in the same manner with a realis- 
tic reaction description. 

Measured cross sections for transfer reactions, such as 



(d,p), have traditionally been analyzed by assuming that 
the reaction occurs in a single step, so that the degree of 
configuration mixing can be gauged by comparing mea- 
sured differential cross sections to those calculated us- 
ing the Distorted Wave Born Approximation (DWBA). 
A spectroscopic factor is extracted, which for (d,p) in- 
dicates the extent to which the final state resembles the 
ground state of the core coupled to a neutron in a particu- 
lar single-particle state. This formalism relies on nuclear 
scattering potentials to describe the interaction between 
the two nuclei in the entrance and exit channels. These 
potentials are traditionally constrained by measurements 
of elastic scattering. In the case of reactions on unsta- 
ble nuclei, the elastic scattering data do not always exist, 
leading to ambiguities in the DWBA analysis. Addition- 
ally, as the deuteron is relatively loosely bound (by only 
2.2 MeV), it may break up in the presence of the target 
nucleus. This deuteron breakup channel can couple to 
the transfer channel, affecting the spectroscopic factors 
extracted in a non-trivial manner. 

To account for this mechanism, Johnson and Soper |20| 
devised the ADiabatic Wave Approximation (ADWA), 
which uses nucleonic potentials and explicitly includes 
deuteron breakup. An extension of this method to in- 
clude finite range effects (ADWA-FR) was introduced by 
Johnson and Tandy [21]. The validity of FR-ADWA has 
been demonstrated in [55], where FR-ADWA was com- 
pared in detail to exact Faddeev calculations. Therein 
it is shown that including the first term of the Weinberg 
expansion works extremely well at these low energies and 
not as well at energies used in fragmentation facilities. 

The two sets of existing ^°Be((i,p)^^Be data, measured 
at Ed = 12 MeV [H| and Ed = 25 MeV [3], have been 
reanalyzed multiple times [TUl dH [23j yielding widely dif- 
fering S values. Even using consistent analyses for both 
sets of data did not result in compatible spectroscopic 
factors [TT1|53|. It is unclear to what extent these discrep- 
ancies in 5* are due to experimental issues (such as the 
normalization of the two sets of data), or insufficiently 
constrained theoretical analysis. 

To help resolve this situation, an extensive series of 
measurements have been made using a primary beam of 
^"Be ions. Elastic and inelastic scattering data were ob- 
tained simultaneously with neutron transfer at four en- 
ergies, covering almost the entire energy range between 
the (d,p) measurements described in the literature. The 
data presented here are for the one-neutron transfer {d,p) 
reaction to bound states and deuteron elastic scattering; 
transfer to resonant states and further scattering data 
will be presented in future publications. The transfer 
data were analyzed using both ADWA-FR [lUlMI, and 
with the traditional DWBA method. 

Data were collected during two separate experimental 
campaigns at the Holifield Radioactive Ion Beam Facil- 
ity [25] at Oak Ridge National Laboratory, using similar 
light-ion detection setups. The initial run was performed 
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FIG. 2. Differential cross sections for transfer to the ground 
state of ^^Be for equivalent deuteron energies of 12.0 (a), 15.0 
(b), 18.0 (c), and 21.4 (d) MeV. Cross sections were calculated 
using the ADWA-FR of Johnson and Tandy [HI El], using the 
CH89 j26 and K-D _27_ optical potentials. Calculated cross 
sections are scaled using the indicated spectroscopic factors. 



at a if'Be beam energy of 107 MeV {Ed = 21.4 MeV), 
with the subsequent runs at 60, 75, and 90 MeV {E^ — 
12, 15, and 18 MeV), with approximately 5 x 10^ particles 
per second. The '^'^Be ions were accelerated from a cesium 
sputter ion source using the 25 MV tandem accelerator. 
Contamination from ^°B was reduced to less than 1% 
by fully stripping the beam ions and tuning the energy- 
analyzing magnet for Z = A. Deuterated polyethylene 
targets with areal densities of 94, 162, and 185 /xg/cm^ 
were used. 

The angles and energies of light-ion ejectiles were 
measured using the Silicon Detector Array (SIDAR 
[H]) (covering 138° < Sjab < 165°), and the first full im- 
plementation of the Oak Ridge Rutgers University Barrel 
Array (ORRUBA [21]) (45° < 6'iab < 135°). The OR- 
RUBA position-sensitive silicon strip detectors (1000 /im 
thick) were mounted at a radius of 76 mm at labora- 
tory angles forward of 95° and 87 mm at more backward 
angles. 

Light ions emitted at forward laboratory angles were 
identified on the basis of their differential stopping power. 
An angular resolution of better than 2° in polar angle was 
achieved. For the purpose of normalization, the prod- 
uct of target areal density and integrated beam exposure 
was determined for the transfer data using the elastically 
scattered deuterons measured in the forward angle OR- 
RUBA detectors, with reference to a low-intensity run 
where the beam particles were counted directly. Protons 
from the {d,p) reaction were detected in the SIDAR ar- 



ray with an energy resolution of w 70 keV at all beam 
energies. The energies of protons emitted from the {d,p) 
reaction at the lowest three beam energies were too small 
to be measured in ORRUBA. However, proton angular 
distributions for both bound states were measured in OR- 
RUBA at Erf = 21.4 MeV with an energy resolution of 
« 200 keV. In the first run, beam particles were counted 
using the new Dual Micro-Channel Plate (DMCP) de- 
tector for heavy recoil detection. A new fast ionization 
chamber, similar to that described in reference |30) . was 
used in the later runs for beam particle counting and 
identification. 

Ground-state angular distributions of protons emitted 
from the ^°Be(d,p)^^Be reaction are compared to ADWA- 
FR predictions normalized to the data in Fig. [2j Opti- 
cal potentials from Varner (CH89) [5B] and Koning and 
Delaroche (K-D) ^27] were used for both the entrance 
and exit channels. No significant differences are found in 
the shapes of the calculated angular distributions. Good 
agreement with experimental data is seen for both of 
the bound states, for all four energies. Additionally, the 
data were compared with DWBA calculations (not shown 
here) , which described the shape of the angular distribu- 
tions well. All transfer calculations in this work were 
performed with FRESCO [31], and adiabatic potentials 
were obtained with a modified version of TWOFNR [32] ■ 
A fixed standard radius and difluseness r — 1.25 fm and 
a = 0.65 fm were used for the bound state. The Reid 
interaction |33j was used to obtain the deuteron wave 
function and in the transfer operator. 

Spectroscopic factors were extracted for each state at 
each beam energy using both the DWBA and ADWA- 
FR formalisms. The results are shown in Figj3} Panel 
a shows the sensitivity to the deuteron optical potential 
(Satchler (Sa) [34] versus Perey and Percy (P-P) [35] . 
keeping the proton potential fixed (K-D) [27]); panel (b) 
shows the sensitivity to the proton potential in the exit 
channel (CH89 [2S] versus K-D with the deuteron poten- 
tial P-P). Spectroscopic factors for the ground (excited) 
state are shown on the left (right). 

The DWBA analysis is sensitive to the choice of op- 
tical potential and there is variation in the value of S 
extracted at each of the four energies using the same 
optical potential. This is most apparent at the highest 
beam energy for the first excited state. These problems 
indicate shortcomings in the DWBA prescription (as dis- 
cussed below). In the case of the ADWA-FR analysis, 
only nucleon potentials are necessary; panel (c) of Fig. 
[3] shows the results obtained with CH89 versus K-D. In 
this case, the sensitivity to the chosen optical potential 
is reduced, and the S extracted for the first excited state 
at the highest beam energy is brought into agreement 
with the results at lower energies. The average S ex- 
tracted from our data are 0.71(5) for the ground state 
and 0.62(4) for the first excited state. 

The inconsistencies arising in the DWBA analysis 
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FIG. 3. (Color online) Spectroscopic factors for ^IBe, extracted from (d,p) data with DWBA (panels (a) and (b)) and ADWA- 
FR (panel (c)) formalisms. Panels a and b illustrate the effect of choosing different potentials in the entrance and exit channels, 
respectively. The error bars are from experimental uncertainties only. See text for further details. 
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FIG. 4. (Color online) Elastic scattering differential cross- 
sections for equivalent deuteron energies of 12.0 (a), 15.0 (b), 
18.0 (c), and 21.4 (d) MeV, including the data of Auton [H] 
and the present study, and optical model calculations using 
the parameters of Fitz [36] , Satchler [34j , and Percy and Percy 

EH. 



come in part from the deuteron optical potentials, as seen 
from elastic scattering. Fig. |4] shows the current elastic 
scattering data compared to those from Auton [8 . It 
should be noted that Auton normalized the data to opti- 
cal model calculations of the deuteron elastic scattering 
using the Sa potential [31] . In the current work, the elas- 
tic scattering data were analyzed using three deuteron 
optical potentials: Fitz [36], Sa [Sj, and P-P [Sg. The 
potentials produce quite different angular distributions, 
and no calculation is able to adequately describe the elas- 
tic scattering data at all four energies. This indicates 



that another mechanism, possibly breakup, is playing 
an important role which the deuteron optical potentials 
fail to describe in the d+^^Be reaction. Furthermore, 
as the Auton data were normalized to one of these cal- 
culations, the spectroscopic factors extracted from those 
data are subject to a significant systematic uncertainty. 
When the Auton deuteron elastic scattering data taken at 
Ebeam = 15 MeV are scaled to the present data (shown as 
stars in panel b of FigJ4]), excellent agreement is reached. 
An alternate analysis of these data using a Continuum 
Discretized Coupled Channel approach has been com- 
pleted and will be presented in a future publication. 

In summary, elastic scattering and neutron-transfer 
measurements, performed with a primary ^"Be beam 
incident on deuterated plastic targets at Ed= 12 to 
21.4 MeV, have been used to study the structure of ^^Be 
halo states in a systematic manner. The poor reliabil- 
ity of deuteron optical potentials for these reactions is 
reflected in the scatter of the spectroscopic factor S ex- 
tracted from the transfer data and its dependence on the 
choice of optical potential. Including deuteron breakup 
through ADWA-FR reduces both of these problems and 
results in a more reliable extraction of S. The average 
S extracted from our data using the ADWA-FR formal- 
ism are 0.71(5) for the ground state and 0.62(4) for the 
first excited state. The quoted uncertainties are solely 
from experimental considerations. The improved energy- 
independence of the extracted S using the present data 
and analysis is in stark contrast to extractions from pre- 
vious transfer data. 
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